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Abstract 

Recent observations of the Galactic center with ground based gamma-ray in- 
struments (CANGAROO-II, VERITAS, H.E.S.S.) have revealed a TeV (lO^^ gy) 
gamma-ray signal consistent with the position of SgrA*. The derived luminos- 
ity of the signal above 1 TeV is a few 10'^^ ergs~^ which is slightly more than 
e.g. the gamma-ray luminosity of the Crab nebula and a plausible identification 
could be a conventional albeit strong gamma-ray source. The observations with 
the H.E.S.S. system of Cherenkov telescopes have decreased the solid angle sub- 
tended by the error box by more than a factor of hundred with respect to the 
previous observations of the VERITAS and CANGAROO groups. After studying 
the observed energy spectrum and angular distribution of the excess as seen by 
the H.E.S.S. experiment, a massive Dark Matter candidate with a minimum mass 
of 12 TeV(90 % c.l.) and an upper limit on the WIMP density for r < 10 pc of 

-I Irx 

1261 Mq pc~^ X (y{(Tv) /3 • 10~^^ cm^s~^) is required to explain the observed flux 
by an annihilation signal. The angular distribution of the excess events is consistent 
with a cuspy profile with p{r) oc with a > 1.0 at a confidence level of 90 %. 
Even though the mass and the cross section of the Dark Matter constituents are 
unexpectedly high in the framework of most models of nonbaryonic Dark Matter, 
it can not be ruled out. 
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1 Introduction 

The presence of a Dark Matter halo of the Galaxy is observationally estab- 
lished from stellar dynamics [1]. At present, the nature of Dark Matter is 
unfortunately unknown but a number of viable candidates have been advo- 
cated within different theoretical frameworks mainly motivated by particle 
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physics (for recent reviews see e.g. [2,3]) including the widely studied models 
of Supersymmetric (SuSy) Dark Matter [4]. On the cosmological scale, the 
presence of non-baryonic and predominantly cold Dark Matter as it had been 
suggested in [5] is well established by many observations most notably of the 
WMAP experiment [6]. 

The implications of the possible presence of SuSy Dark Matter in the Galactic 
center has been studied in numerous publications including predictions for the 
7-ray flux expected from annihilation of weakly interacting massive particles 
- WIMPs (see e.g. [7,8,9,10,11]). With the advent of the next generation of 
Cherenkov telescopes a gamma-ray signal from the direction of the Galac- 
tic center has been established: The recent claims for gamma-rays from the 
Galactic center by the CANGAROO collaboration with a signal S = 9.8 a 
[12] and a possible excess claimed by the VERITAS collaboration on the level 
of 3.7 a [13] have already stimulated first speculations on the interpretation 
of the signal as annihilation radiation from Dark Matter [14]. Finally, the re- 
cently commissioned H.E.S.S. system of Cherenkov telescopes has detected a 
signal from the vicinity of Sgr A*[15]. The different confidence regions of the 
detections by the different instruments in comparison with a 21 cm radio map 
[16] are indicated in Fig. 1 and the flux measurements are given in Fig. 2. 
Note, the EGRET measurements taken from [17,18] give the range of possible 
flux values for non-resolved objects in the gamma-ray background dominated 
region and should be interpreted as upper limits for a contribution of Sgr A*. 
The upper limit indicated in Fig. 2 at 1 GeV is from a dedicated analysis of 
the high energy EGRET data [19] assuming a photon index of 1.5. 

At first glance, the differences in the reported TeV flux by the CANGAROO 
and H.E.S.S. groups might suggest variability at the lower end of the energy 
spectrum thus ruling out an exclusive annihilation origin and hinting at a 
more conventional compact gamma-ray source related to Sgr A*. However, 
given the positional uncertainty of the CANGAROO II and VERITAS results 
evident in Fig. 1, it remains possible that in fact more than one (possibly even 
variable) source may exist in the direction of the Galactic center. Interestingly, 
the observations of H.E.S.S. indicate the presence of even further sources in 
the field of view (see Fig. 1 of [15]). 

Further observations with H.E.S.S. and CANGAROO HI in the coming years 
will help to clarify the situation which might be consistent with multiple 
sources which are not resolvable by single telescope observations. Instrumen- 
tally, the CANGAROO and H.E.S.S. observations have shown discrepancies 
in the flux values for other steady sources: SN1006, PSR B1706-44, Vela, and 
RXJ1713-304 [20,21,22] - a similar situation might be present in this case. 

The VERITAS energy flux given in Fig. 2 is converted from the 40 % Crab flux 
as measured with the Whipple telescope, assuming a photon index of 2.2, and 
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using the Crab flux as determined by H.E.S.S. to compare the two instriiments' 
observations. The error bars are given as the sum of statistical and systematic 
uncertainties. The two measurements are in reasonable agreement and given 
the possible positive bias in the flux estimate for a source detected at the level 
of less than 4 standard deviations should be considered to be consistent within 
the uncertainties. 

The discussion here is focussing on the results obtained with the H.E.S.S. sys- 
tem of telescopes which achieves a superior angular resolution in comparison 
to the previously pubhshed results as is evident from Fig. 1 ^ . The recon- 
structed source location is, within the errors, consistent with the position of 
Sgr A*(space angle = (18 ± 42)") and appears to be consistent with a 
pointsource for the given angular resolution of 5.8' (50 % containment radius 
[15]). 

The H.E.S.S. observations (spectrum and angular distribution) arc interpreted 
in the framework of a Dark Matter annihilation scenario to investigate the 
consequences for the WIMP mass, density, and spatial structure and to un- 
derstand whether such an interpretation is feasible. 



2 WIMP mass 



In order to quantify the minimum mass required to explain the observed energy 
spectrum in terms of an annihilation signal, an assumption on the shape of 
the annihilation continuum spectrum is required. In the very generic case of 
a self-annihilating massive WIMP, it is assumed that the branching ratio for 
qq, W+W-, and is quite large. 

Subsequent hadronization leads to the production of charged and neutral 
mesons which decay finally into electrons, neutrinos, and gammas. The gamma 
signal is described by a fairly hard differential energy spectrum with an ex- 
ponential cut-off d^/dE oc E"^-^ ex-p{—E/Ec) with Ec ~ O.lm^ [7]. However, 
the exact shape of the spectrum depends on details of the branching ratio of 
other processes. Namely, for xx ~^ 't^i direct radiative losses produce a second 
hard component. 

For the purpose of calculating a lower limit on it is assumed that the 
branching ratio for xX ^ or M^^M^^is close to one. The minimum mass 
requirement is somewhat relaxed if xX ~^ '^^ is of importance. 



These data were taken during the build-up phase and will eventually be super- 
seded by the performance of the full telescope array 
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Fig. 1. Position of the TeV emitting region overlaid on a 21 cm radio map [16] tracing 
mostly non-thermal (synchrotron) emission features. The bright spot inside the 
H.E.S.S. confidence region (marked as a white circle) is Sgr A*whereas the extended 
ring like feature to the east of it is Sgr A*East. Note the difference in the accuracy 
of the different instruments. 

Fitting a generic 7r° decay spectrum to the data results in a 90 % c.l. lower 
limit on the mass to be larger than 12 TeV. 

In most models that attempt to explain non-baryonic Dark Matter in the 
Universe, the relic particles are WIMPs produced thermally during the early 
phase of cosmological expansion. A very reasonable candidate for this particle 
is provided by SuSy models which naturally predict the existence of a stable (if 
R-parity is conserved) Majorana particle that can self-annihilate. Generically, 
this SuSy particle is referred to as the lightest stable particle (LSP). 

The detection of gamma-rays with energies up to 10 TeV as seen with the 
H.E.S.S. telescopes implies that the WIMP mass has to be sufficienly large 
to explain the observed spectrum. The surprisingly large mass inferred from 
the observations appears at first glance to be inconsistent with cosmology and 
the general conception that the mass of LSP in SuSy models should be close 
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to the electro- weak unification scale [23]. However, it has been argued that 
a multi-TcV WIMP is in fact natural in Minimal Supergravity models [24]. 
The upper bound on the possible mass of a WIMP that decoupled thermally 
during the early phase of expansion is around 40 TeV [25] for a cosmological 
non-baryonic Dark Matter density of O^dmh^ = 0.1. 

The average annihilation cross section {av) for SuSy WIMPs is in the right 
range to explain today's Dark Matter density as determined by the WMAP 
mission: ^dm.h'^ = 0.1 16 ±0.020 [6]. However, given that for the thermal decou- 
pling Qdmh^ oc {crv)~^ and that (av) oc "rn^^, rather tight constraints on the 
mass of the WIMP < 600 GeV can be derived [26]. However, as already 
pointed out in [26] there are specific manifestations of SuSy models (even in 
the constrained, universal models) that allow for larger WIMP masses which 
still produce the right relic density: 

(1) For specific SuSy models in which the LSP is almost degenerate in mass 
with other SuSy particle (e.g. f) coannihilation depletes the relic density 
even for a small annihilation cross section. 

(2) For a non-thermal decoupling process, the relic density as it is today is 
not depending on the cross section. 

(3) Relaxing the universality of the Higgs multiplett masses at the GUT scale: 
In this (and in the general SuSy) case, it is conceivable that a massive 
(multi-TeV) WIMP with the right relic density may exist. Thess WIMPs 
are generally almost pure Higgsino type Neutralinos with a large cross 
section. 

In the first case, the cross section for annihilation remains small and therefore 
requires high relic density values in the Galactic center to explain the observed 
7-ray flux. In the special case of a degeneracy between the LSP and other 
SuSy particles, careful tuning is required to achieve the right relic density. 
For a massive LSP beyond a few TeV, even co-annihilation does not suffice to 
reduce the relic density. A non-thermal scenario invoking Q-balls [27] which 
replenish the Dark Matter content via evaporation at some time after the 
decoupling has been studied in the literature [28,29]. In this and other non- 
thermal scenarios the cross section for WIMPs could be rather high. 

The last case is more attractive but again requires tuning of the Gaugino 
masses at the unification scale which might appear less attractive from the 
point of view of a natural theory. 

Generally, a massive WIMP as suggested by the observations with a suffi- 
ciently high cross section is not the obvious candidate in the framework of 
MSUGRA but could be the LSP in slightly less general theories. In this way, 
thermal decoupling would produce the right relic density and the observed 
flux could be consistent with a Dark Matter density which is not violating 
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other astrophysical bounds. 



Ultimately, a neutralino mass exceeding 40 TeV is excluded by the unitarity 
bound in the framework of thermal decoupling [25] . 



3 WIMP density 



Given the measured flux, it is possible (without assuming a specific halo shape) 



to derive an upper limit on the product {pdm)r<Roy (c^^) This is truly an upper 
limit, as can be seen from the following argument: 

Assuming that the true density distribution is spherically symmetric and fol- 
lows a general form [7] 



^^''^ (r/a)^[l + (r/a)^](/^-^)/" ^'^^ 

The halo density profiles predicted from N-body simulations of hierarchical 

structure formation follow a power law roughly between 7 = 1 [30] and 7 = 1.5 
[31]. The annihilation flux observed from within the Galaxy is proportional to 
the line of sight integral 



J(^) = JoX J dlp^ (^■^R^ -21 R cos + (2) 

l.o.s. 

averaged over the solid angle observed: 



{J{0))AnAn = Jo X JdQ Jdlp^ (3) 
with the commonly used normalization for the Galactic Center: 



8.5kpc(0.3 GeVcm ■^)^ 



cm^ (4) 



For reasons of simplification and in order to calculate an upper limit, it is 
convenient to use a truncated version of Eq. 1 which drops to zero outside a 
spherical volume with radius Rq. Obviously, the derived density to match an 
observed flux would be higher than in the case of a continuous halo distribu- 
tion. For a truncated density function with Rq much smaller than the distance 
of the observer to the sun, the line of sight integral can be safely replaced by 
a volume integral and the observed differential gamma-ray flux is given by: 
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d$ I I dN (av) I 2 
- X dV 



dE 2 47rrf2 dE ml 



^{p'W (5) 



Considering a truncated density function which follows a single power law up 
to a radius Rq dropping to zero for larger radii: p{r) — po x {Rg/r)" with 
a = 0...1.5. 



JdVp 



3 — a 



and for (p^ 



6 — za 



which gives for the ratio (p^)/(p) 



2. 



(P^) _ (3 - a? ... 
(p)2 3(3 - 2a) ^ ^ 

This ratio is > 1 for a > 0. For o; = the trivial case of a constant density 
and (p2) = (p)2 is recovered. For a > (p^) > (p)^ which implies for a 
given flux (oc (p^), see Eq. 5) a lower average density than in the case of 
a constant WIMP density ^ . Therefore, calculating an average density with 
a — from the observed flux sets an upper limit on the actual average density 
for a realistic radial density distribution. 

For convenience, Eq. 5 can be rewritten with appropriate units: 



'^=0.54xl0--(cm^s)- ^^'^ 



dE ' ' ' (5000GeVc-2cm-3)2 3 x IQ-^ecm^s-i 

fm^Y\^ (3-a)2 ^ / d \\ f Rq VdN 



\TeYj ^ 3(3 - 2a) ^ [s.5 kpc J ^10 pc J dE 



2 For the sake of brevity, this is only shown for < a < 1.5, but it holds true for 
any a > 0. 
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In the more general case of a radial density profile extending beyond Rq, the 
line-of-sight integral of Eq. 3 will naturally give larger observed fluxes (oc 
( J(0))An Ail) than for a truncated density therefore relaxing the requirement 
of the average central density (p) even more. In this sense, the approach of 
calculating an average, truncated density is truly an upper limit on the average 
density as expected for a realistic, smooth Dark Matter halo distribution. 

The parameterization for the decay spectrum is of course depending on the 
nature of the WIMP and cannot be given in general. For the case of a very 
massive WIMP considered here, the dominant annihilation channel is into a 
W~^W~, Z^Z^, and qq pairs. A parameterization of the form 



is a reasonable description of the decay spectrum (with E and in TeV, 
£;o = 2-10-4TeV). 

In order to calculate the upper limit on (p) it is assumed that the entire ob- 
served flux is produced by annihilation of WIMPs. The H.E.S.S. data-points 
are taken from [15] with the assumed emission region being covered by the 
point-spread function of ~ 5.8'. The actual limit on the source size is only 
3'. However, this assumes a Gaussian surface emission proflle, whereas a vol- 
ume source with constant cmissivity is expected. The upper limit on such an 
emission proflle will be somewhat bigger than 3' but smaller than the point 
spread function. The upper limit depends on the assumed emission region as 
ip) oc Ro'^\ 

The H.E.S.S. data points (open triangles) are indicated in Fig. 2 together 
with the CANGAROO (open boxes), VERITAS (open diamond), and different 
EGRET flux values (open symbols from [17], filled symbols from [18], and 
upper limit from [19]). A best fit model of the H.E.S.S. data requires 



m;^ = (18.9 ±3.7) TeV 

GeV / (av) \ "^^^ a/3 (3 - 2a) 



= (1261 ± 124) ^ X ( fc^ 

^ ^ pc3 \^3 X 10-26cm3s-V 3- a 

with xfedid.o.f.)^lA3{13). 

The confidence regions given by the minimization for the two parameters 
(p) and is given in Fig. 3. The best fit position is indicated by a small 
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cross. 



For the sake of completeness, a similar procedure can be applied to the CAN- 
GAROO data giving the following values: 



m;^=(l.l±0.3) TeV 



(P)r<47pc = (7840 ± 1170) ^—r X 



GeV / M /3(3_2«) 



c^cm^ V 3 X 10 cm^ s ^ / 3 — a 



- (206 ± 31) X ( <'^"> V"' V="^-'"' 

with Xredi^-^-f-) ~ 0.81(4). The value given for the density is higher than 
calculated in [12] because of the factor of 1/2 in Eqn. 5 that has been neglected 
in the past [32] . As pointed out in the Introduction, the CANGAROO-II result 
is not consistent with the H.E.S.S. observations hinting at either a mismatch in 
the cross-calibration of the two instruments or a variable object which might 
coincide with one of the more conventional gamma-ray sources close to the 
Galactic center. 



The required average density depends on the cuspiness of the halo profile 

and scales according to Eqn. 9 with ^^3(3 — 2a/ (3 — a) (for a < 1.5) and is 
reduced by a factor of 2 with respect to a constant density halo distribution 
when a ~ 1.4. 



4 Spatial WIMP distribution 



Given the excellent angular resolution of the H.E.S.S. telescope arrays of 5.8' 
for individual events [15], spatial extension of the source down to this level 
can be probed. The signal appears to be point-like (i.e. of spatial extension 
smaller than the point spread function) with an upper limit of 3' for a source 
with a Gaussian surface brightness profile. 

The expected source extension from Dark Matter annihilation depends cru- 
cially on the shape of the radial density profile. For a quantitative study, a 
simplification of the radial density profile given in Eqn. 1 is applied: 



p(r)=po X 



{r < a — 8.5kpc) 



: (r > a) 



(11) 

fl2) 
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Fig. 2. A summary of data and best-fit models for WIMP annihilation from the 
Galactic center: H.E.S.S. (open triangles), CANGAROO (open boxes), EGRET 
(solid and open circles), 10m Whipple telescope of the VERITAS collaboration 
(solid diamond). 

Since the observations cover only the inner 0.5° of the Galactic center (equiv- 
alent to 75 pc), the actual value of a is not crucial. 

The observed surface brightness is calculated by folding J{ip) (Eqn. 2) with the 
average point spread function of H.E.S.S. for the observations approximated 
by a double Gaussian 



fpsfW = fo ■ (exp(-0.5V^V^^i) + 1/8.7 • exp(-0.5V^Va2')) (13) 

with (7i = 0.052°, (72 = 0.136° (J. Hinton, priv. communication), and fo as an 
arbitrary normalization: 



—oo 

with Pq chosen such that / di/jp{ilj) = 1. The resulting angular distribution is 
shown in Fig. 4 for different values of a. Using the values for various a, a 
90 % c.l. (one sided) for a > 1 is calculated which is suggestive of a halo with 
a halo steeper than predicted by Navarro, Prenk, and White [30]. 
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Fig. 3. For the fit of the model energy spectrum to the data confidence regions 
are indicated for different confidence values. The best-fit model is indicated with a 
small cross. 

5 Discussion 



TeV gamma-ray emission from the Galactic center has been observed by dif- 
ferent different ground-based instruments. The CANGAROO group claims a 
signal at a significance close to 10 a and the VERITAS group sees a indications 
for a signal observed in archival data. Finally, the H.E.S.S. system reported 
an overall excess of 11.8 a from the Galactic center region from observations 
during the build-up phase. 

However, taking the H.E.S.S. data alone, no indication for variability is seen 
and moreover, the source location is consistent with the position of Sgr A*(A^^ = 
(18 ±42)" with a systematic error of ~ 20") and an upper limit for a Gaussian 
surface density source of 3'. 

Taking these results, a hypothetical Dark Matter annihilation scenario is dis- 
cussed. The observed energy spectrum is reasonably well described by the 
continuum emission from SuSy WIMPS (neutralinos) with a mass of 18 TeV 
(lower limit of 12 TeV) putting it far out of reach for future accelerator ex- 
periments. 
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Fig. 4. For different values of the slope of the density profile, the expected observed 
angular distribution (after folding with the point spread funetion) allows to separate 
the case of a < 1.0 and a > 1.0. From the good agreement of data with a point 
source, a steep (cuspy) halo is required. 

The required average Dark Matter density is ( 1092 ± 107) -MqPc^^ (for (av) ~ 
3.10-26 cni^s"^ and a — 1) for the inner 10 pc of the Galaxy with a cuspy profile 
p oc with q; > 1 in order to fit the observed angular distribution. The upper 
limit on nonbaryonic Dark Matter density is close to the stellar (baryonic) 
mass density which is p*(r) ^ (1900 ± 700) x (r/10pc)-i-^±°-i MqPc-^[33]. 
The total inferred mass M(r < 10 pc) = (3.1 ± 0.3) x 10^ Mq is well below 
the stellar mass of (1.5 ± 0.6) x 10^ M© in the inner 10 pc. The rather high 
central values of the density is consistent with the formation of a mini- cusp as 
a consequence of interaction with the central dense stellar environment [34,35] 
and baryonic compression [36]. 

The synchrotron emission emitted by the electrons/positrons injected through 
the annihilation process constrains the maximum slope of the radial density 
distribution [37,38] which rules out a spike close to the black hole [39]. The 
constraint is further relaxed if fast diffusion or convection removes the elec- 
trons/positrons from the vicinity of the black hole. 

The angular resolution in gamma-rays achieved with H.E.S.S. allows to probe 
the actual shape of the Dark Matter profile for the Galactic center. With more 
data taken with the full telescope array operational, a further improvement on 
the constraint on the cuspiness is expected and test the different possibilities 
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for the origin of the signal. 



At this point, the considered annihilation scenario of WIMPS as the origin 
of the observed flux of gamma-rays is not ruled out. The other detectable 
channels of annihilation i.e. excess of positrons and anti-protons in Cosmic- 
rays are not affected by a mini-cusp scenario as outlined above. The large 
mass of the WIMP may also produce a signal for neutrino-telescopes but the 
sensitivity of the current instruments is not sufficient to detect neutrinos from 
the Galactic center [40]. In principle, bounds from big-bang nucleosynthesis 
could be of importance but as outlined in [41], a very massive WIMP would 
not have a strong impact on abundances of light elements. Future observations 
of variability in the gamma-ray flux or a possible extension would be of great 
importance for the interpretation of the signal. 

The author is grateful for discussions with F. Aharonian, J. Buckley, A. Ja- 
cholkowska, J. Primack, P. UUio, H.J. Volk, and members of the H.E.S.S. 
working group on Dark Matter. 
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